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Analysis of transverse impact behaviour of 
simply supported and built-in beam is investigated in 
the present work, with and without taking into account 
the contact deformation. The impact is assumed to be 
purely elastic. Hertz's law for contact deformation 
is assumed to hold, 

Experiments have been conducted to verify the 
theoretical prediction. Non-dimensional impact force, 
and maximum strain have to be plotted against impact 
parameter which depends on geometries and the elastic 
properties of the beam and striker. 



CHAFTBR 1 


INTRODUCTION 

1.1 TRANSVERSE L.PACT ON BEiii'': 

During the last four oecades, scientists aiic! 
engineers have exhibited an increasing interest in the 
solution of problems concpiuied with the dynamic loading 
on structures, particularly in the case of impact. 

Actual examples of this type of loading may be found in 
the fields of tool design, foundry and machine shop 
operations, protective ordnance, vehicle accidents and 
in many other areas. 

But the severe mathematical complexity encoun- 
tered in treating the problem as a whole restricts the 
study of impact phenomena to impact involving only 
simple type of geometry^ and in the case of machine or 
building structures the simnlest and most important type 
of structural element would be a simply-supported or 
built-in beam subjected to a transverse impact load. 

The duration of impact, the maximum impact 
force exerted and peak values of stress or deflection 
are the quantities of interest for such a situation. 

In the neighbourhood of impact point the transient 
stress distribution is obviously three dimensional, 
the waves are not only propagated along the beam but 
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are also reflected from all surfaces. At more distant 
stations the stress will be more uniformly distributed 
across the width, exhibiting the well known bending 
features with stress variation at a given cross section 
occuring primarily as a function of depth. This pattern 
also extends to the region of contact point after remo- 
val of load. Because of analytical and experimental 
limitations the three dimensional theory does not 
appear to be feasible. 

In order to make the problem amenablo to 
theoretical treatment it is assumed that the variation 
of stress across the width can be neglected. Thus the 
problem boils down to a two dimensional problem of plane 
stress. Additional simplification is brought in by 
neglecting the rotary inertia and transverse shear .. 

1.2 historical DSVRLOPMEHT 

In the field of one dimensional transverse 
impact on beams many scientific research workers have 
contributed various analytical and experimental approa- 
ches. Various mathematical models to the problem have 
been suggested since 1800 for tli'^ solution. The simple 
approximation for elastic deflection of a beam under 
transverse horizontal impact obtained by equating the 
energy of striking body to the resilience of beam may be 
attributed to Young |lj. Hodgkinson £2]| analysed the 
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problem with simple consideration of momentum transfer 
and justified them with experiments. Ills results were 
given a theoretical basis by an important analytical 
investigation made by Cox |J3 j to explain the results 
of a Royal Commission's exneriments [4j • Later a 
complete solution was presented by Saint Venant and 
A.Flament l_5j who recognised the presence of vibratory 
motion which materially influenced the form of elastic 
curve. The nineteenth century investigators attempted 
to rationalise experimental results with vibration 
theories by postulating the existence of a quasi-elastic 
layer near the contact point* This concept was super- 
ceeded by the theory of local contact deformations 
introduced by Hertz £6^ which has found wide use in 
spite of the static elastic nature of its derivation. 

Realizing that the local deformation of the 
contact region was an important influence, Timeshehko [7 J 
devised a solution combining the contact pressure equa- 
tion of Hertz with the lateral vibration theory. Solving 
the differential equation of forced vibration under a 
central force varying arbitrarily with time, he related 
the beam deflection and local de'f'ormation to the striking 
masp displacement* His extensive work paved the way for 
further investigations with more modifications in the 
mathematical models and verifications by the experiments. 
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Zener £8^ calculated the energy losses during impact 
and gave an experimental verification. A further simpli- 
fication in the analysis for the impact problem was done 
by Lee in his Doctoral Thesis. He assimed 

that because the duration of impact is small compared 
to the fundamental time period of the beam, first mode 
of vibration of the beam is sufficient to be considered 
for the analysis. But in some cases, particularly for 
impacts with large energies, the results obtained did 
not confirm the earlier theory given by Timoshenko. 

Mason l^lOj showed in his experiments the influence of 
multiple collisions by the striker, particularly v/hen 
the striker mass or the striking velocity is increased. 

In the field of photostress analysis extensive work was 
done by Foppljjlll to verify the result with Timoshenko's 
Theory. Galerkin's technique and Collocation' method 
have been used by Eringen [12') in deriving apnroximate 
and quickly solvable results. 

During the last decade experiments were conduc- 
ted by Horst Schwieger j^lD , l4J with steel spheres 
impacting on beams. He assumed that because the impact 
duration is much less than the natural time period of 
the beam, in the vicinity of contact point the pealis 
of stress and force occur much earlier than the arrival 
of the reflected waves and so they are unaffected by 
the boundary conditions. 



1.3 PRESENT WORK 
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In thp present investigation two different 
independent mathematical methods have been suggested 
for the analysis of the impact problem. The first one 
does not consider the contact deformations while the 
second one takes into account the relative approach of 
the bodies during the process of contact. 

Experiments have been performed to verify the 
tijeoretical predictions. Effects of impact velocities, 
the mass of the striker, boundary conditions etc. have 
been studied for various beams. Applicability of both 
the mathematical methods ha s been studied comparing the 
results predicted by them with the experimental values. 
Non-dimensional quantities such as maximum impact force, 
maximum stress have been plotted against a non-dimensional 
parameter (A) in order to predict quickly the maximum 
stress developed in a structural element subjected to 
an impact with a l-oiown mass and striker velocity* 

Influence of various parameters of beams and the striker 
(such as the striker velocity and geometry) on the value 
of impact parameter X has been studied. 



CHAPTER 2 


ANALYTICAL INVESTIGATION 

2.1 mathematical models 

In the analytical investigation t ¥0 models are 
considered. In the first the effect of the contact 
deformations of the structural element due to the 
striker mass is neglected. This implies that the 
displacement of the striker and the team at the impact 
point are assumed to he the same (fig.l, model l) , 
that is, 

Wq (|, t) = ¥ (t) (2.1.1) 

¥here Wq = Displacement of beam at impact point 
¥ = Displacement of striker mass 

i = C = striking point on the beam. 

In the second model the effect of local contact 
deformations of the beam due to the striker is taken 
into account, that is, there is a relative motion bet¥een 
the striker and the beam during the impact. This 
relative motion (fig.l, model 2), ¥ " is given by 
the Hertzes contact deforuiation equation as follows 
for purely elastic impact. 

3/2 

F(t) = K [w(t)«W^(| 5 t)J (2.1.2) 

where K, the Hertz constant y depends upon the material 
properties of the two bodies and their geometries* For 



the Case of cylindrical shaped striker X'fith lieLiisoherlcal 
end and the beam, it would be [”15^ 

K = I ® ^ r2 (2.1.3) 

where = Eg = B = kodulus of e la:3ticit37' for bean 

and striker in ate rial; 

= ’^>2=1)= Poisson’s ratio for both the 
bean and the hanner and 
r = radius of the hemispherical end of the 
hamiuer. 

The following assumptions have been made in 
both the mathematical laodels s 

(a) Natural frequency of free vibra.tions of the 
striker is very high* Therefore its any type of 
effect is neglected* 

(b) Collision is perfectly elastic tliat is no perma- 
nent Indentations occur in the two colliding 
bodies at the point of contact* 

(c) The supports are ideal* They do not absorb any 
part of energy during the impact* 

(d) Presence of only one collision during the jjmpact 
is assumed. For higher initial velocities, 
multiple collisions talie place and they affect the 
occur ranee of peak values of impact force or 
strain (fig. 2). 


s 


2.2 ANALYSIS V-flTrf NO CONTACT DYFOR-iiiTION 

An approximate method is developed here 
wherein the beam vibration equation is solved aiid the 


applied impulse is replaced by a su.its'hle initia.l 
condition j 5, 15J , The equation of motion for the 
transverse vj.bration of a beam js 

S I + fA = 0 (2.2.1) 

o x^ d t'^ 

where = Young's modulus of elasticity for beam 

I = Moment of inertia of the beam cross section 
'f - Mass density 
A = Cross section of the beam. 

The solution of the above equation can be 


written as 


ffO 

■h i 


W(x,t)= / A.Sin^x +B.Cos^x +C.3in h Vx 


i=l 


Cos h I X 


Ej_Sin <A[_t+F^Cos 

( 2 . 2 , 2 ) 


where 


^ = Eigen frequency, and 


A _ SL- 
^ - fA 

Substituting the following boundary conditions 


for simply supported end conditions in eq. (2,2,2), 
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W(0, 

t) = 0 


(2.2,3) 





0 ¥ 
a x '2 

(o,t) = 

0 

(2.2.4) 

h ¥ 

<5 X 

('2,t) = 

0 

(2.2.5) 

The constants 

Bi = D. 

= 0 and, 



Ci = -A^Cos(-l|^;)/Cos h(-S^) 

With the initial condition 


the value of = 0. And finally tjie 6q.(2.2.2) reduces 
to the following forrn 


W(x,t)= 


where, G 



(2.2.7) 


( 2 , 2 . 8 ) 


and, 

X^=Sin( J^x)/Cos(-=^)- Sinh(5’^x)/Cosh(~i-) ( 2 . 2 . 9 ) 


Substitution of eq.(2,2,7) in the following boundary 
condition 








( 2 . 2 . 10 ) 
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yields the characteristic equation in 0^ as follows, 


0^(tan 0^--tanh 0^) = 2 M 
where M = Mass ratio = M^/m^ 

i = 1,2, 




b.L 


cxD 


( 2 . 2 . 11 ) 

(2.2.12) 

(2.2.13) 


Evaluation of of the eq,(2.2.7) im^olves 
replacenent of the impulsive force by an initial velocity 
condition for the beara. It \rill be assuxned here that 
the striker imparts a velocity indistinguishable from 
Vq to an infinitesimal section of the beam just at the 
point of contact. Then, momentum of this section and 
striker should be equal to the initial momentum 
that is, 




(2.2.14) 


where integration is carried out with respect to the 
total mass of the beam and striker. 

From eq. (2,2,7) , 


i. I. T. KANPUR, 

i. , 




<’>0 


|-|(x,o)= 

1=1 


Am* 

„ ^ <».« Mis'll 


(2.2.15) 


or 


oo 


^ Gj_ X^ = Y (2c) 

i=l j = 1,2,, 


(2.2.16) 


oO 
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If eq,(2.2,i6) is multiplied hj m 2 and 
L 

evaluated at x = 2 and added to the in'ce^ral of 

eq.(2.2.i6) with respect to beam mass dm-.=m. /L dx, 

1 

there is obtained in view of symmetry, 

- L/2 


f 2m. 


Z s 

^ J 

O 

1=^1 


X, X, an: + 10 X (^) X.(-) > 


i D 


0 2 


L/2 
2m_ r 


J 

o 


X. t|'/-(x) dx + mg X. (|) (2.2.17) 


which is 3^st the representation of integral in 
eq,(2.2.i4). Utilising the orthogonality of X^ and 


X. , 




L/2 

L 

L o 


J X^ y'(x)dx + nig X^Cg) y^(L/2) 


(2.2.18) 


2m. 

T 


L/2 

f o 2 

X. dx + mo X, (L/2) 
i 2 1 


o 


^ w 

'h t 

'h w 


Applying the initial conditions at t = 0, 

(2,2.19) 


= y(x) 


for X / •- 
2 


= Vq = ^(L/2) for X = ~ 


eq,(2,2.l8) is reduced to the following form : 


( 2 . 2 . 20 ) 



^2 Vq (L/2) 


2m 2 2 

r X^ dx + m^^ (L/2) 
0 


4V 


o 


0 


1 


JL 


2M 


(2.2.21) 


^Lcos^ 0^ Cosh^ 0j_J 


Finally substitution of eq.(2,2.2l) in 
eq.(2.2,7) gives for simply supported end conditions, 




W(x,t) 


^2 

L 


a 


0 . 

i=l ^ 


20^x 20.x 

Sin — Sinh — ~ 


.L ^ 


Cos0J“ Cosh 0. 


■SH 


Sin 


.,/y 2 2 

40^ a 


^^0320^ CoS?0i 02 


( 2 . 2 . 22 ) 


Similarly for clamped end conditions the 

expression for deflection is 
o oO 


V_L 

W(x,t) = ^ 


i=l 


0 , 


2^ . 


( l-Cos0^Cosh0i) (Cos0^-Cosh0j_) (Sin0^+S±nh0j_) 
_(l-Cos0^Gosh0.)^+ M(Cos0^-GosH0.)^ 


1 


uJ 


Sinh 


20^x 


it 

' mm 


Sin 


20.x 


Cosh 0j[_« Cos 0j^ 
•) 


4^ a^t. 


20,.x 20.x 

Cosh -Cos 

Sinh0^+ Sin0i 


Sin( 


(2.2.23) 
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where 0^'s are determined from the following charac- 
teristic equation, 

(l-Cos0.Cosh 0.) 

M = -i 3: (2.2.24) 

Sin 0j_Cosh0^+Cos0j_Sinh0j_ 


Strain is obtained by differeirciating eqs.(2,2.22) and 
(2,2.23) twice with respect to z and then substituting 


it in following equation 


(: (x,t) 


H ^ ^ 

2 ^)x2 


(2,2,25) 


The impact force F(t) is evaluated by knowing 
acceleration from eqn.(2.2,22) and (2,2,23) and then 
multiplying it by mass of striker mg as follows 


F(t) 



(2.2.26) 


Numerical methods and results are discussed 
in Chapter 4, Unfortunately convergence in the case 
of eqs,(2,2.25) and (2,2,26) is very poor and results 
have been provided with consideration of first ten 
harmonics. 


2,3 ANALYSIS OTH CONTACT DEFORMATION 


In this section the local contact deformations 
of the beam and the strilccr are taken into account. The 
response of the beam and the impact force are calculated 
using 


i) 

ii) 


the Hertz’s Law, 


e q. (2.1,2) 


the Dynamics of Hammer 


1 

J 
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From the available standard methods [l5, iSj 
for the solution of beam subjected to a force varying 
arbitrarily with time and acting at distance x = C, the 
deflection of beam is given as follows s 

cO . 


W(x ,t) 


X, 


q Cos(u) t) + Sin(vO t) 
o i i i 


i=l 


where 


Xi(C) 



o 


t 

j F( .Sinq) .(t- 't^dX 

j IL 

o 

(2.3.1) 


u) ^ = Eigen frequency of beam 

X. = Eigen function of beam 
1 

J = Mass density of beam material 
A = Uniform cross-section of beam 


(2.3.2) 

(2.3.3) 

(2.3.4) 

(2.3.5) 


oO 

q 7 X. .= Initial displacement of beam 
° i ' 

oO 

X./oO.= Initial velocity of beam 

i=l 


(2.3.6) 


(2.3,7) 


For the beam initially at rest, expressions 
(2.3,6) and (2.3.7) become zero, and eq*(2.3.i) 


reduces to 



OO 
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W(x,t) = 


X,^(C) 


F (X* ) Sin t- xO d X' 


i=l 1 ° 


(2.3,8) 


For a simply supnorted beam, 


i j2 


iif! (SL_)1/2 

O V •» r. n_ 


'djh 


(2.3.9) 


« Sin (ifTx/L) 


(2,3.10) 


X^'^dx = L/2 


(2 *3 .ll) 


and force acting at x = C = L/2, the eq.(2,3,6) yields, 


iSln^ (|^) 


i=l ^ 


tr 

J FCxOSin%(t- tOd-X’ 

Imm 


(2,3.12) 


Rearranging the above equation following one 


is obtained * 


Wo(i-« 


u oc.> 

r r 

F ( X ')-1 > 

ifMt ^ 


X Slnul^(t-V)la 

CjOj 

(2.3.13) 


The value of Sin^(i IT /2) fluctuates between 
the limits +1 and -1. Thus representation of summa- 
tion sign in eq«(2.3,l3) in the form of integration is 





possible [ 1^3 or, 


-ro = i j ®iii-i±i20 di 


(2,3,14) 


where i means continuously varj^ing quantity of tlie 
states of vibration arid is given by 

i = L/( I /2) = u>L/('TT.C ) (2,3,15) 


where 

= phase velocity 
= cO(i) 

* v” 

From eqs.(2.3,l5) and (2,3,9) 


. d<jO (■dfhs^'^ du> «L fdfh')!/^ o'^^d 


to Ti 


(2.3,16) 


The integral in eq,(2,3,i4) becomes 

f Sjn . 1 . (^ 13 ) 1 / 4 . a oj 

J Ld -n Bi 

o 

I- /d'fh^l/^ / ,xl /2 

= :p^*(^“) ./2 7 t , (t-i:*) 

Substitution of eq, (2.3.17) in (2,3,i3) gives 

= ji] ^ • 


(2.3.17; 


■yi SI 
1 /2 

( t-tn dX ^ 


_ /2_ j 

“V^ • djh • ^ El 


1/2 

F(q::*).(t- 't') d 

(2,3.18 
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Defining a beam parameter 


I ^ * - 1 


^ SI ^ 


eq.(2,3,i8) becomes 
t 

Wo(|,t) f F( d'c> 


J 

o 


From the dynamics of hammer, 


d^w 

Dig = - F (t) 


dt‘ 


vrith the following initial conditions ; 


w 


= 0 


fc=o 


(2,3.19; 


(2,3.20: 


(2,3.21) 


(2,3.22) 


dwi 


= V, 


t=o 


(2.3,23) 


Integrating twice the eq.(2,3,2i) following equation is 
obtained for the displacement of hammer 


w(t) = V t 
o 


t 


r 


S 

o 


F(X*)-(t-T*)d Z' 


(2.3.24) 


The two eqs.(2,3,20) and (2.3.24) are -elated by 
0 q.( 2 ,i, 2 ) as 

3 /2 

F(t) - K “* ^0 (2,3.25) 

where K = Hertz constant defined in eq.(2,l,3). 
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Subs'bi'b'u.'ti.jig eqs«(2,c','20) ario. (2,1 
eq.(2.3.25) 


.24) in 


-2/3 2/3 


F = v^t - r F( r').(t- r') . 


o { 
✓ 
o 


1 


L 


d 


Defining 

2/5 


;j275~“'~][75 “ Reference Time 
K ,1^ 


i 0<* ^2 

= "iT7S” 


= Impact parameter 


and using the dimensionless variables 


'T = 

^ T 


tf = FwCt) - VK|,t) 1 /T 


follovring integral equation Is obtained, 

j-_r- 

0 -x 


(2,0.26) 


(2,3,27) 


(2,3,28) 


(2.3,29) 


(2.3.30) 


<^=X- \ r) d - - ;\ l ^)i/2 ^ 


(2,3.31) 


Proceeding on similar lines equation can be 


obtained for curvature at x = L/2 as follows : 


L 






(2.3.32) 
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and final expression for strain will bo 






u 


O 

(2,3,33) 

The quality in thr:' integral sign in eq.(2,3,33) is 
denoted by f ( 't , A ) i.e . 


f(ir,A) = 


3/2, . 
.5“ (hr- 


- 1/2 _ 
X) d ''d 


(2.3.34) 


.j 

o 

Squation (2.3,31) is solved numerically by 
finite Difference Method [jL7~|j laiowing Impact force 
is calculated from eqs.(2,3,30) and (2,3,25) as follows: 


F = K(T (2.3.35) 

Deflection strain 9(^ ,t) are obtained 

respectively from eqs.(2.3.24) and (2,3,33) and (2,3.34) 
by the numerical procedure of integration a'^’ter knowing 
0^ or F. 



CHAPTER 3 

EXPERI14SHTAL ARRMC-EMEHT 


3.1 BRIEF DESCRIPTION OF THE SET UP 

An experiment has been des-lgnec^ to verify the 
theoretical predictions. The exrerimontal arrangement 
is given in Figure (3,1 Details of the set up are sho^m 
in plates 1 to 4. Beams of uniform rectangular cross 
sections have been tested. Table 1 gives the list of 
tests conducted with various hammers and beam cross 
sections. End supports are shown in Figures (4,5)1 
they have been so designed that the loss of energy 
during the impact is minimum. This is accomplished by 
avoiding any sort of deforinations at the support. To 
achieve this hard steel spheres and hardened steel flats 
have been used in the supports. 

The hammer is supported by means of low stiffness 
stainless steel tube to minimise the energy absorption 
and also to ensure the movement of the striking mass in 
a vertical plane. At the top, the suspension tiock 
[fig.[6jlhas been precisely made to avoid frictional 
forces. Sliding arrangement at the top makes possible 
the adjustment of the suspension point of the hammer 
wherever desired. Scale attached to the top supporting 
block measures (of the order 0,5®) the angle of swing 
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for the haEimer, Neglecting the friction losses and 
assuming pendulum to be a simple one, expression for 
striking velocity is obtained by equating the ma/Cimun 
potential and kinetic energies of the ha-Uior j 

1/p 

\ = [2gR (1 “ Cos e)] C3.1.1) 

where R = Length of Pendulum. 

3 . 2 INSTRUI^NTAT ION 

In studying experimentally the response of the 
beam subjected to transverse impact, various parameters 
such as bending strain, deflection of the beam and the 
acceleration of the har/mier during the impact have been 
measured. Fig. (7) shows the line diagram for the perma- 
nent record of outputs in the storage oscilloscopes and 
their simultaneous triggering with the acceleration 
signal. The very short duration of the impact (of the 
order 100 micro seconds) needs effective design of the 
electronic circuitry as well as the proper selection of 
the transducers for measuring various parameters 
involved. In the following sections transducers together 
with their electronics have been discussed, 

3.3 ACCELERATION MEASUREMENT 

For shock measurements following two types of 
accelerometers have been used because of their linear 
behaviour for a large range of acceleration and other 
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reasons discussed later. For Jetalled specification 
see Appendix , 

1. Endevco accelerometer. Model 22i5 
Working range 0 to 2000 g 

Voltage sensitivity = 14,55 pk mV/pkg 

2. Endevco accelerometer. Model 2225 
Working range 100 g to 20,000 g 
Voltage sensitivity ,50 pk mV/pkg. 

These accelerometers using ■oiezoelectric 
crystals have fairly high natural frequency of vibra- 
tion (of the order 20,000 cps) and therefore they are 
suitable for very high response. The circuitry involved 
is shown in Figure (8) . 


Following is the formula used for calculating 
the voltage sensitivity of the accelerometer (^iSl s 


E = Ecal 


4- c, 


^p ^cable’*' ^amp 


(3.3.1) 


where , 

®cal ~ Voltage sensitivity of the accelerometer 
with standard capacitance 
Cp = Accelerometer capacitance 
Ccabii Cable capacitance 

^amp “ Amplifier capacitance including the 


residual capacitance-. 
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3.4 STRAIN i'iEASUREl'iSNT 

It is obvious that the maximum bending strain 
occurs at the middle point of the beam when the load 
acts at that point. Therefore as sliov.n in Figure (3) 
two gages A and B have been mounted on the back side of 
the beam at the middle point. They are connected to 
the Wheatstone bridge as shown in Figure (9). Calibra- 
tion of gages is done by shunting fixed resistances 
across the active gage |~3 8 . The equivalent strain 

is given by 

n ^ .i (3.4.1) 

equivalent Fg ^ 

where 

Rg = Gage Resistance = 120.6 + .2 JTu 

Fq. = Gage factor = 2.06 + .02 
6 

R = Calibration resistance 
c 

For three different calibration resistances 

following are the equivalent strains using eq. (3.4.1) 

i- =28 micro strain 

■ 1 

R^=10®_n. 

f =56 micro strain 

^ 2 

R^= .51X10® 

c 

'3 

R^=.2xl0® xm. 


l43mlcro strain 



24 


These values are plotted against the correspon- 
ding voltage output of the circuit in Fig. 1. And the 
strain sensitivity is found to he 

= 116,80 pk micro strain/pk mV 

Temperature compensation is provided using the 
dummy gages in the other two arms of the bridge. 

3,5 D1IFLECTI0N nEASUREMEITT 

The deflection will occur at the middle 

point of the beam for the case under consideration. So 
it is of main interest to measure this value. Solar 
cell, which utilizes a beam of light to measure the 
displacement is used because it has 

(a) "fery high response because of the absence 
of any kind of mechanical inertia effects. 

(b) high resolution without any noi^e in the 
output signal on the oscilloscope screen. 

(c) linear behaviour of the solar cell for a 
fairly good range of displacement. 

The circuit used for the solar cell is shown 
in Figure (lO). The calibration is shown in Figure (12) 
and the sensitivity of the transducer is found to be 

Eg = 0.00042 in/mV 

For specifications of the transducer see the 


Appendix 



CHAPTER 4 


RESULTS AI® DISCUSSIONS 
4.1 theoretical RESULTS 

L 

The deflection of the oeam at z = 2 
obtained from eqs.(2.2,22) and (2.2.23) using IBM 7044 
digital computer. The convergence of the series was 
very good and the first ten harmonics were sufficient. 
These results have been plotted in Figures (l5, 19, 23, 
27, 3i) for all Cases and it is obvious that maximum 
deflection obtained for clamped case is always smaller 
than that for simply supported case. But unfortunately 
the high sensitivity of second derivatives with respect 
to X or t of eqs. (2.2,22) and (2,2.23) makes it nece- 
ssary to include more harmonics in the calculation. 
Particularly in case of beam -with clamped ends, conver- 
gence in case of force or strain calculation is very 
poor 1^12, lO’J , Inclusion of higher harmonics involves 

greater computational efforts in handling very large 

l20 

quantities (of the order e )• Therefore forces and 
strains were calculated only for simply supported beam 
taking into account only the first ten harmonics. 

For the second analytical approach (Sec, 2 .3) 
finite difference technique was used to solve for the 
non-dimensional impact force ( ry ) and non-dimensional 
Strain (f) from eqs. (2,3,3i), (2,3.34). Knowing C^and f 
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other quantities like % (|, t), (- (J, t) , F(t) and 
duration of impact can he obtained. Results for all 
cases (see Table l) have been shown in tables 2 to 6* 

In Figures (l3 to 36) all the theoretical results (force, 
strain, displacement, time duration) together with 
experimental results have been plotted against 
or time. 

4.2 EXFERIMEINTAL RESULTS 

For the five cases mentioned in Table 1, 
acceleration of hammer or the impact force ^ accele- 

ration) were recorded for different striking velocities. 
The results have been given in Figures (13 to 36) ajid 
tables (2 to 6), 

The value of B for each beam has been calculated 
knowing the n-icural frequency of fundamental mode of 
vibration for simply supported end conditions of the 
beam : 

^2 1/2 

= — [si/dPhJ (4,2,1) 

Table 1 shows o'.’s and the corresponding E for 
each beam, 

4.3 GRAPHICAL PLOTS 

Figures 11 and 12 are the calibration curves 
for strain gages and solar cell. For each case of the 
beams (Table l) four plots corresponding to maximum 
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impact force, mazimum strain, mazimiini deflection and 

impact duration have been drawn against the striki-Qg 

velocity V^CFigs. l3 to 32), The non-dimensional 

, 3/2, 

maximum impact force ( ) and maicimuii: strain 

max 

(fjijg^) are plotted against impact parameter ( -^ ) in 
Figures 33 and 34. To give an idea of the difference 
in the magnitude of forces and strains and their beha- 
viour with respect to time for the two iiiatheniatical 
models 1 and 2, plots (figs. 36 and 26) are made for 
a particular case ^t. of the striker = 2,92 lb ; 

= 3.5 in/sec 5 beam used = beam 3 aJ« 

4.4 EFFECT OF IMPACT VELOCITY 


Increment in impact velocity Vq increases the 
Various quantities like impact force, strain, deflec- 
tion as evident in Figures 11 to 32. The isjjact dura- 
tion time has a tendency to become lower with increasing 
velocities as shown in Figures 16, 20, 24, 28, 32, Since 
in the mathematical model d escribed in section 2,2 no 
energy is dissipated in the contact deformations, it is 
expected that force, deflection and strain should be 
greater as compared to those given by analysis in 
section 2,3 or the experiments. This is seen to be 
true from Figures l3 to 32 except for the case consi- 
dered In Figures l3 and 17 which corresponds to Beam 1 
and Beam 2 respectively. Thus it can be concluded that 
beam dimensions do effect the impact force aS calculated 
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from eq,(2,2.26) for simply supported end conditions. 

And it can "be said that for higher values of depth and 
height of the beam the results obtained from eq.(2.2.26) 
for the impact force are better. In the case of Figures 
16, 20, 24, 28, 32 only the experimental values and 
those obtained from eq,(2,3.35) are plotted for comparison 
because it is obvious that in the first model the impact 
duration will be much less as shown in Figures 35, 36. 

4.5 EFFECT OF BOUNDARY CONDITIONS 

The effect of the boundary conditions on the 
various parameters was studied for the cases of Beam 1 
and Beam 2. In Figures 13 and 17 the impact force 
obtained from eq.(2.2.26) is lower then the experimental 
values for simply supported and clamped boundaries* 

In the case of strain (figs. 14, 18) the experimental • 
values for simply supported ends are very close to those 
corresponding to eq.(2.2.25) but for clamped ends, the 
experimental values are lower than those for simply 
supported beam showing the effect of boundary condition. 
Similarly in case of deflection the efioct of boundaries 
is clear from Figures 15 and 19. In all the above cases 
the experimental values for s.1mply supported beam are 
20 to 40 percent higher than those for the clamped beam. 
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4.6 EFFECT OF PLASTIC DEFORMATION 

This Gffsct is found to b© considerable in 
Case of iripact force. At low impact energies the expe- 
rimental values of F^^^ are in the nen ghbourhood of the 
predicted ones, but in case of higher impact energies 
the experimental values show a remarkable deviation 
from the predicted values as seen in Figures 13, 17, 

21, 25, 29. In the beginning with lower values of Vq, 
experimental curve has a tendency to follow the curve 
obtained from eq.(2.3.35) which takes into account the 
contact deformations ; here force obtained from eq*(2,2.26) 
need not be compared because of the obvious deviationsj, 
but later for higher values of the slope of experi- 
mental curve reduces. This decrement in the slope is 
mild in Beam 1, Beam 2 and Beam 3 a but appreciable in 
Beam 3B and Beam 3C^ This deviation may be attributed 
to the local plastic deformation. This shows that more 
accurate theoretical analysis is needed in the case of 
large impact velocities. 

4.7 EFFECT OF IMPACT PARAiffiTER A 

Figures 33 and 34 shov; the dimensionless quanti- 
ties (maximum impact force) and f , „ (maximum 

strain) plotted against the impact parameter A respec- 
tively. For a particular configuration of beam and 
hammer, the value of X depends upon the velocity Vq* 
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Both the predicted values and the experimental values ha"' 
been compared. In case of Figure 04, the agreemeat is 
fairly good, but in Figure 33, at some stations the 
deviations are considerable. This may le attributed to 
the following reasons ; 

(a) The force indenta.tion relation as given 
by eq,(2,i,2) is valid only for static 
compression of two isotropic elastic bodies. 
In dynamic case it may change . 

(b) The relation (2,1.2) is valid for smooth 
curved surfaces. Tven a little amount 
of inaccuracies involved in the fabrica- 
tion of the hammer may change the value 
of K considerably. 

(c) The supports cannot be ideal. 



CH AFTER 5 


CONCLUSIOKS Am SUGGSSTIOIIS 
FOR PURTiiER WORK 


5.1 CONCLUSION 

The theory discussed in section 2,2 does not 
take into account the contact deformations and therefore 
it Can be treated as more idealistic as compared to the 
other which considers the effects of contact deforma- 
tions. The results obtained from it were also on the 
higher side as compared to those obtained from experi- 
ments or the other theory. Moreover it needed more 
computational efforts in calculating the second order 
derivatives with respect to time or x-coordinate • The 
trouble of convergence creeped in when the series was 
truncated at the fifth harmonic. If more than ten 
harmonics are considered the trouble of handling large 
quantities is encountered. Ultimately computation was 
done with only ten harmonics# In the case of simply 
supported beams the convergence was satisfactory, but 
for the clamped beam the situation was worse. However 
deflections in both the cases were obtained without 
any trouble. 

The second theory (sec. 2,3) gives better 
results as well as tbe computational effort is also 
less. Finite difference technique for solving the 
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non-linear integral equation worked very satisfactorily. 
Moreover because of the dimensionless quantities involved 
in the computation, the .job was simplified to a greater 
extent. Thus it can be concluded that the second theory 
which takes into account the contact deforniations should 
be preferred for use in case of the transverse impact 
on beams inspite of the assumption that the impact is 
purely clastic. 

For high velocities of impact a I’emarkab.le 
departure was obtained in case of the predicted results 
of impact force from the experimental values. This may 
be accounted to the limitau on of the theory that nature 
of imnact should be purely elastic. In the absence of 
this requirement the Hertz theory needs modifications 
in the calculation of the value of K, because it may 
not remain constant and may be influenced by the aiiiount 
of indentations which occur in the beam. It is seen 
from the experimental results that the boundary condi- 
tions effect the beam resuonse. 

Instrumentation involved in the measurements 
of various parameters needs very careful attention. In 
the Case of acceleration measurement, number of trials 
with Various similar accelerometers were made to check 
their calibration and reliance for very high response. 

The response of the gage largely controlled by its 
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inertia, was assumed to be sufficiently high to permit 
the recording of the strains during impact the duration 
of which was about 200 micro seconds ^ 20 J . Using the 
solar cell for the measurement of deflection proved to 
be a very good and precise method in displacement measu- 
rement and also suitable for very high resoonse because 
of the absence of any typo of mechanical inertia effects. 

5.2 SUGGESTIONS FOR FURTHER WORK 


Regarding the analysis, improvement is possible 
if one takes into account the variation in the Hertz 
constant and the effect of permanent indentations in the 
two bodies. For exaaple the force indentation relation 
can be expressed as 


n 


r 


(5.2.1) 


F(t) = K [wCt) - Wq t) 
where n in present work was 1.5. 

Here the solutions may be investigated with 
the value of n other than 1.5. Instrumentation should 


be improved for very high impacts particularly in case 
of more rigid beams where the rise time of stresses or 
forces reduces to a very low value. Also the effect of 
multiple collisions should be studied. Attempts should 
be made to simplify the analysis to predict better 


re suit s 
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T/iBLE 1 

LENGTH OF BEiXIS = 40.125” 


No. 

5 D 

5 in. 

'Size 
f H 
JL„in. ^ 

S rad/ 
sec„^ 

!l P H S 

5 ^5 psi 

5lb/in'^5x 10® 

5 rap 

$ 

5 lb 

ilHaranier 

5 radius 

5 In. 

Beam 1 

1.985 

♦483 

174 

.282 30.0 

2.92 

2 

Beam 2 

1.957 

.956 

130.5 

.271 29,8 

2.92 

2 

Beam 3 a 

1.995 

.976 

362 

.283 30.5 

2.92 

2 

Beam 3B 

1.99S 

.976 

352 

,283 30,5 

3.76 

2 

Beam 3C 

1.996 

.976 

352 

#283 30 #5 

8.95 

3 
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Beam 1. ; J' =.282 lb/ in®, D = 1.985 in, 

H = .483 in,, E = 30xl0%si, nig = 2.92 lb, 

= 10.85 lb, M = 3.72, OC= 3.3911, lb“^in. see 
K = .3l051xl0®lb in.~^'^^,v= .3 


in/ sec. 


‘TxlO’ 

sec. 


li 


Sk&Yq)' 


5 max 


‘FExp. 


7 

96.711 

2,606 

.412 

547 

.623 

.522 

.750 

14 

84.194 

2,792 

.396 

1255 

.691 

,460 

.646 

17 .,6 

80.521 

2.855 

,391 

1640 

.582 

.412 

.628 

21 

77.640 

2.908 

.386 

2045 

.574 

.351 

.660 


28 73.304 2.993 


.380 


2880 


.561 


376 


615 
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T/^LB 3 

Beam 2. DrJlA ; .282 Ib/in.^, D = 1.957 in., 

fi = .356 in, B = 29.8 x 10® psi, m 2 = 2.92 lb., 
™1 ~ lb., M = 2.59,o(=: 3.4453 lb“^in.sec ”®^^5 

K = .308 X 10® Ib.in"^-^^, V = .3 



3.6 110,300 4,030 ,3l5 240 ,441 .400 .600 

7 96.070 4.318 .300 550 .416 ,436 .466 

10.6 88,600 4.497 ,292 896 ,402 ,322 .416 

14 83.640 4.630 .286 1266 .392 .285 .434 

17.6 79.990 4,730 .282 1655 .385 ,267 .407 
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t;j5LS 4 



Beam 3.1. 

H = .976 

DiH;. ; ,283 

in., B = 30x10 

lb. /in.® 
^ psi,®i[* 

, B ~ 1 • 

.013387 : 

995 in. 

Ib^^in, 

5 

-3/2 

sec 


'l}' = .3, 

m = 22. 

2 lb, m 

2 = 2.92 

H 

m 

li 

7.6, 



K - .3157x10® lb 





Iji/sa 

-6 

5Txl0 ^ 


(jf 3/2 
■'max. 




3.5 

110.360 

,840 

.709 

240 

1.222 

.846 

1*480 

7 

96.070 

,899 

. ^690 

651 

1.184 

.590 

1.380 

10.6 

88,690 

.,937 

,680 

895 

1.162 

.458 

1.362 

14 

83.640 

,964 

,673 

1265 

1.147 

.363 

1,342 

17.6 

79.990 

.986 

.667 

1654 

I#l35 

.320 

1.355 
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5 

Beam 3B . DATA : 3 = .283 D = 1.995 in., 

^ S 3 

H = .976 in., B = 30x10 psi,c>(= .013387 l'b~^in.sec. 

= .3, m = 22.2 lb., mg = 3,76 lb., M = 5.9, 

K = .3157 X 10^ lb. in.“®'^ 
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Ti'^Ln 6 

Beam 3C. DATA s P= ,283 Ib./ln.^, D = 1.995 in., 

H = .976 in,, E s 30x10® psi, 6( = «0l3387 Ib’^in.sec^'^^ 
.= .3, ='22.2 lb., nig = 8.95 lb., M = 2.48, 

K = .387 X 10^ Ib.in."^'^^ 



3.5 

159,300 

2,140 

.461 

7 

138,660 

2.295 

,443 

10,6 

127.860 

2,389 

.433 

14 

120.720 

2,459 

,425 

17.6 

115,460 

2.515 

.420 


610 

.717 

.871 

1.398 

1171 

.684 

.504 

.923 

1900 

.664 

.388 

.821 

2690 

.650 

.306 

.793 

3510 

.639 

.263 

.720 
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APPENDIX 


1. SPECIFICATIONS OF TRANSDUCERS 

(a) solar cell 

Size 1.1/8'' X 1.1/8” x 3/l6” 
Capacity, Volts .3- to .45 V 
Capacity, liA lO-lo 
Cat. No. Sill 

(b) ACCELSROI'ffilTSR 

i) Endevco, Model 22i5 

Voltage sensitivity with 300 pF 

external capacitance =14*5 mV/g 
Capacitance 10850 pF. 
ii) Endevco, Model 2225 

Voltage sensitivity with 300 pF 

external capacitance =.5 mV/g. 

(c) STRAIN GAGES 

Make - Rohits Co, (India) 

Internal Resistance 120,6+' •! ohms 
Gage factor 2,06 + .02 
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2, FINITE DIFFERENCE J^fflTHOD FOR SOLVING EQ.(2.3.3l) 

Eq,(2.3.3i) is a non-linear Volterra integral 
the general form of which can be given as flT'J 

X 

y(x) = g(x) +J'K|x,t 5 yCt)^ dt (1) 

o 

On setting x = 0, h, 2h ....... *rh in succe- 
ssion in eq.(l) and using the trapezoidal rule for 
numerical integration the value of y(x) at station 
is given as 

^r+1 " ^r+l ^r+l ^r+l’ ^r+1^ 

( 2 ) 


where 




K (x D X 5 Y ) I 
r+l’ s’ s J 

s=l 

Eq,(2) is solved iteratively by means of the 
procedure ; 

y(p+ 1)=: g _ + U + -I h K(x 5 9 Y ^ ) (3) 

r+l °r+l r+l r+l’ r+l r+l 

(p) th 

where p = 0, 1, 2, and Y^^^ P appro- 

ximation to Y^^^. 

The following Predictor-Corrector method with 
only one iteration process gives good approximation to 
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the value of Y with an error of the order of (h )* 
r+1 


Predictor 


"^r+l ®r+l ^r+1 ® 


r+l° "r’ ^ 


Corrector 

o 

Y = a -I- U + hK (x ^^3- X„J., 3 
r+1 °r+l r+1 r+l"^ r+1’ r+x 

( 5 ) 


where has been defined in eq. (2)* 
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COMPDTBR PRO G RMMEi-LISTIHG 


cc characteristics roots for sihply supported case 

OliMEfsiSlON PHIS(SO) ,PHIC(50) 

A fv] = E e A 8 

PYE=A„^(-ATAN( lo ) 

X1S=1«5 

X23=ioA 

Xl = 2o 
X2=2ol 

N = 40 

DO 120 1 = 1 >N 

121 Th1=SINH(X1S)/C0SH(X1S) 

Th2=SINH (X2S) /COSH(X2S) 

F0 = 2,*AM-XlS^MTAN(XlS)nTHl) 

FN = 2 „ *aM-X2S^'- ( T an ( X23 ) -TK2 ) 

Xf\l = FN*XlS/ ( FN-FO ) +F0-X'X2S/ { FO-FN ) 

X1S=X2S 

X2S=XN 

IF ( AB3{XN-XlS ) «LT o «001 )GO TO 122 
GO TO 121 

122 PHIS(I)=XN 
X1S=PHIS ( n+PYE 
X2S=XlS"^®l 
PUNCH 888, PHIS U) 

120 CONTINUE 

CC CHARACTERISTICS ROOTS FOR CLAMPED BEAM 
D0123I=1,N 

12A CZl = lo-COS(Xi)-«-COSH(Xl) 

CZ2=i«-COS(X2)*COSH(X2) 

FO = AM-Xl*CZl/ (SIN(X1 ) *C0SH ( XI ) +GOS ( Xl ) *S I NH I X 1 ) ) 
FN=aM-X2*C22/ ( S in ( X2 ) ^f-COSH { X2 ) +COS ( X2 ) *SINH 1X2}) 

X N = F N X 1 / I F N — F 0 ) + F D * X 2 / ( F 0“ F N ) 

X1=X2 
X2 = XN 

I F { ABS ( XN-Xl ) oLT» »OOl)GO T0^23 
GO TO 12A 

123 PHIC(I)=XN 

XI = PHIC( I )+PY£ 

X2=Xl+0.1 
PUNCH888 ,PHI C( I ) 

123 CONTINUE 

saa FORMiAT ( lX,El3.>b ) 

STOP 

END 

sentry 


A4 
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CC 

cc 


12Q 

L21 

888 


5 


i^itTSiOu Fur SpLViNu tuSs'lisi 

A r'i = i-i A -b S RATIO *• i“i i / M 2 


bcU* ( 2 , 2 5 


uiFiLRSlOi^ PHIS( Ay ) »PHIC<AO) sFSibO) >FC( 80 } 

U I rA:-'.ibION UEFS ( 3 & ) > DEF C ( 50 ) j SIRENS ( 50 ) sSTREMC { 50 ) 

DIMENSION DEGREE (6) 

data DEGREE/1, ,2o ,3o ,5, ,6./ 

AM=2oA8 

SM = 8oy5/(32e2^«-12o) 

D0120I=1 ,A0 


READ 888 ,PHIS( I ) 


001211=1,40 


READ888,PHIC( I ) 
FORMAT! IX, ElPoP ) 
PYE=4o*ATAN{ 1. ) 


AL=40»125 

Do 6 1 J= 1 ,6 
D05l=l,50 
SIRENS! I )=0. 

ST-ENC! I )=0» 

DEFS!I)=0. 

FS(I)=0® 

FC! I )=n- 

DEFC(I)=0o 

TW=0. 

THETA=PYE*DEGReE( U) / 180 . 

V0=SQRT{ 2.^32, 2^12. *104, 25*! 1. -COS ! THETA ) ) ) 


PRINT 11 

PRINT 55, VO 

ROW= o283/(32»2*1<i« ' 

0=1.995 


H=.976 

£I=30o5*10«*M6*D4H4*3/12. 

A= ! E I / ! ROW*DMH ) ) ** o 2 5 
PRINT ll 

PRINT 22, A, AM, El , VO 
PR+NT 11 
PRINT 33 
D01I=2,50 
I 1 = 1-1 

FFS=FS! I-i )*16.*A**2*V0/AL**2*SM 
FFC=FC! I-l) *8 .*A442*V0/AL*42 *SM 
WS = DEFS( I-1)MAE^-*2>V0/A4*2 
WC=DEFC( I-l )MAL**2*V0/ !A**2)*1./4 o*2. 
SS=STRENS !T-1 )*2 .*V0*H/A**2 

SC=STRENC ! I-l )4V0*H/A**2 
T = TW*AL*^^2/ (4,*A4*2 ) 

PRINT44, I I ,TW,W3,WC,S6,SC»T ,FFS,FFC 
TW=TW+.0005 
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uO-s:j = i,io 
PHS=PHIS(J) 

PHC=PHiC<J) 

TMT=TAN ( PHS ) -TANH (PHS ) 

TPi = TAN( PHS)+TAlMH(PH3) 

CCP = l*/C03{P + S)*'>4-£;— I,/C0SHvPilo)-x-#z + ii .*Ai''i/PHS** - 
31 = 5 IN ( TW*PhS**Z ) 

^2 = SIN { TA'^i-PHC*-»2 ) 

COCC=COSH{PHC)-COS(PHC) 

5 1 5 1 =S I NH ( PHC ) +5 I iNi ( PHC ) 

iJEFS { I ) = L)tFS ( I ) + TiViT*Sl / ( PHS*-"- J*CCP ) 

^TREN5( I ) =bTK£N5( I }+TPT*Sl/ I PHS*CCP ) 

F S ( I ) = F 5 ( I ) + T i'"'i T -* S i -* P ri S / C C P 

oRi^ = ( i , -CUP ( PhC ) *Cu5H ( PmC i ) -* ( COS ( PriC ) -COpH (PhC) )-*{PIN{ PHC > +S I NH ( P 
IMC ))/((! • — COo ( PhC ) -J^-COSh ( PhC ) ) **2+Ai'i'"- ( LOo C PHC ) —COSH ( PHC ) ) -**2 ) 

PFFC( I )=o£FC{ I )-fdRt\-*p2*( (PI,'>ill(PrlC)— Pli'j(PHC) ) / COCO~(i_Oph(PHC) 
i -COS ( PhC ) ) / SISI ) /PHC*-*2 

•pTRENCI I )=STKENC( I )+bRK*{ ( S I NH ( PHO ) -fP I N I PHC ) ) /COCO- ( COSH ( PHC ) + 
lCOS{PHC) )/SISl )*S2 

FCU )=FC( I )+dR:<*S2-«PHC**2-;i- ( ( S I rviH < PHC J -S 1 N ( PHC ) >/ COCO- ( COSH ( PHC ) - 
iCOS(PHC ) ) /SI SI ) 

2 CONTINUE 

C G0)MTU-^L For OVERFLOW 

D014K=11 »2S 
PHS=PHIS(K) 

PHC=PHIC(K) 

S 1 =S I N ( T W*PHS**2 ) 

S2 = SIN C TW-x-PHC*-^/ ) 

OEFS( I ) =DEFSOI ) -H TAiN (PHS )*Pi/ ( i® / (COp (PHp } } **2+2 .-^mm/PHS-*-k-2 ) /PHP**: 
COCO=COSH(PHC)-COS(PHC) 

Sl5l=5lNH{PHC )+SIN(PHC ) 

PTKEnS( I )=PTRE!'')p( I ) + TAN(PhS)*Si / (1 «/(Oijp(PHS) )**2 + 2, *AH/ PhS**2 ) / P I 
IS 

FS ( I ) =FS ( I ) +TAN (PHS )*pi*PHS/ ( 1 ./ (COS (PhS) ) **l+2 . *AN/PHS**2 ) 
PRK=SInH ( PHC ) *COS (PHC ) / (COP (PHC ) **2+AH ) 

OEFC( I )=0EFC( I )+bRis.*o2*( {SIi\li{PhC)-bIts(PHC) )/COCU— (CL*SH(PHC)— COS 
KPHC) )/SISI )/PHC*-»2 

PTRENC ( 1 ) =pTPvEimC ( I )+dRK* USIrm ( PHC ) +bi N (PHC ) ) /CoCo-(COSH ( PF 0 + 
ICOP(PHC) )/SISI )*S2 

FC( n=FC( I )+bRA*P2-*PHC*-*2-» ( ( S I I'iH ( PHC ) -P IN ( PHC ) )/COCO-(COSH(PHC) - 
1 COS (PHC) ) /SISI ) 
lA CONTINUE 
i CONTINUE 

6 CONTINUE 

11 FORMAT {/iX» 120 ( lH-> ) 

22 Format ( lX»*A=*»Ej.6.b »3X»*AM=*»E16.P ,bX,*EI=*»£iP.3 >3X »*vO=*!.Elo.d 

1 ) 

3 3 F uRi'iAT ( 2X » * I No* » pX » *T iw* > lAX » -**■ wS* > IAa » *i/'sC* » XAX » *pp* » lAX » *SC » 1 + X » 

1*T*) 

AA F0RMAT(2X» IA,SEi3.3 ) 
bb Format ( IX »* vu=* » Elo . 3 ) 

STOP 

End 


SENTRY 
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CC tZK-HLRTZ CONS'lAiNiT jl-ZcFo TiiFt 

, ALpiD= Impact PaRAMCTER jALFA = btAM PARAl-iZTtR 
Sj E X £ C U I E W A 7 F 0 R > R E W I N D 

ilBFTc MAIN 

'0 i hens I ON S I GMA ( S 1 ) , TOW! 5 1 ) ? FTl ( 5 1 ) ? D£F (SI) 

Dimension ptosi ) ,w( si ) 

Dimension vo(6) 

Dimension degreeis) 

DATAOEGREE. / 

R=10A.12S 

Rl=3* 

E = 3 0 » 3 -jf i 0 a ■)<-* 6 
EZK=.666*E*SQRT (R1 )/„91 
0=1.995 
H=.97 d 

ROW = , 283/ ( 32,2*12. ) 

S M = 8 , 9 5 

S M = SM / {32,2*12. ) 

PYE=4,*ATAN( 1.) 

EI=0*H**3*E/ 12. 

DO 3 I = i , 6 

ThETA=PYE*DEGREE { I ) / i80« 

3 VO(I )=SQRT(2. *32. 2*12. *K*( l.-COS{ THETAI ) ) 

DO 19IJ=1,6 

c t-reference time 

c almd-impact parameter 

C ToW-T /REF. TIME 

G Sigma — {vj“WO)/t*vo 

C ALFA — beam parameter 

T= (SM**.4 ) / ( EZR*^ »4*V0 ( IJ ) **.2 ) 

PRINT354T,EZK 

AlFA= ( D*ROVv'*H/E I ) ** . 2 5 / ( SORT ( 2 ,*PYE ) *D*ROW*H ) 

ALiMD= ( ALFA*SM)/SQRT (T ) 

PRINT44,ALMD!>ALFA,90( U) 

TOW(1)=0. 

DO i00l=l>30 
100 TOW{ I+1)=T0W( I )+,1 

SIGMA(i)=0« 

AH = .1 

c -Time interval for tow 

DO 101=2,51 
il=l-i 

C— PREDICTOR 

SUM=0» 

9 SUM=SUM+ ( SIGMA { J ) **1 , 3 ) * ( TOvii t I )-TOVU J 3 +ALMD*' ( TOvv ( I ) -TOw ( J ) ) **,5 3 
CALC = AH*t ,5*SlGlMA(i)**l«5*(TO'W£ I ) -TOW ( 1 ) +ALM0* ( TOWN I )-TOW( 13 3*^»3 


1)+SUM3 

bIGMA( l) = TOw{ I ) -CALC-. 3* AH*S I GMA Cl 13**1 


.5*C TOW C 1 3-TOw( I I 5+ALMD 


ik- 


1 (TOWC I )-TOW( I I 3 )**.5 3 
IF (SIGMAC I 3 .6E.0. )G0 TO 10 
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n r. r, 
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I F ( 5 I (>;'■■ A ( I ) o L T o 0 . ) i'i = I - i 
GO TO 7 

C ti U T H L R El IT IS i'M 0 T N t. £ D c D 

iO C U r\i T I l AU t 
C — n — C/RRECTOR 

C S T IR A 1 1'* A S F R o 0 F T I I'-i t 

7 FTL(i)=0o 
OFF ( 1 ) =0o 
GOBI ,N 

I F ( (5i6RA( I )-SlU:''!rt( 1-1 ) ) abtoOo );-lMX=I 

I I = i - i 

S U i’-'i ~ U o . 

U i'-'j il/' “ o 

D02uJ = l^ II . . ^ u 

oURiD = bUiMO + SIGr-iA(J i Uv.O I ) - Tuw t J I ) - ^ . 3 

^0 SUf'i = SUM+SIGMA{ J ) ( 70W ( I ) -Tuw ( J ) 1 » J ^ ^ 

F T L ( I ) =AH^^ { o I GHA ( 1 ) 1 o 1 OW ( I ) -Twyi. ( i ) ) ® 

UEF( I )=aH*( obvtSIGiMMl 1 ) To.'. ( I ) -Towil ) ) ^ 

8 CONTINUE 

PRINT 11 
UO BO I = i 9 N 
o I i’-i A = S I Gf'lA ( I D ** 1 e ^ 

3U PRINT^2» I sSImAjFTLC n 9T0W( I) 

S I MAX = S I OiRA ( I'iAX ) >H5-i o S 
PR I ivT 3 3 sFlAX 9 S I MAX jFT L ( lArtX ) 


l.S.-r, KANPUR, I 

4,»T« I 





calculation of force .strains UI SPLACbi'itN I 


AO 
b b 

11 
19 
11 
22 
3 3 




OJ 


PRINT 77 
UO AO I = 1 5 i''! 

SIi'iA-SlGi'iAl I * '•■''■■'■“P . 

PT ( I )=EZR-SIi»iA* ( T*VO ( I.J) ^ 

oTl=io73L*tr<Ov</L)-:;-*„b-x-(MLFA/^o)''t_.c.,s vovu; 

ThV|E = TOO( I )-:J-T , .r - u.--v c r t i 

»v ( I ) = LZR* ALF A’>Pv 0 ( 10 ) 1 o J'''‘T '' '"^3 " Ut . 

PRIilT 65 5 I s P ! ( n s S r 1 s Vi ( I ) s T iFiE 

CONTINUE 

faORiMATl iXs IosAE1d„ 3 ) - lO X , > ' , T I N * , ax » ’VOEFL EG T I ON */ » 


v3Xs«-FTL( I )-x,iOXs*REF TIi^iE-"") 


FGR.-.AT ( / xXs*STEP 
CONTINUE . _ , , , 

FOPi lAT ( SX .---sTEP^'- s iAX 5 ■'- ( o I Oi’>A J i o - 

I o D » / s iX s 1 20 ( iH— ) ) ^ t- / 'A = * i-M 6 a '■) I 

FORhaT ( / / sBOX s*T = *s Eioc. j s-a-EZA s 

STOP 

END 

SENTRY 
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Thesis 
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of traa^ers! 
impact on teaJHS. 






